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Abstract. We present some results on the existence and nonexistence of cen- 
ters for polynomial first order ordinary differential equations with complex co- 
efficients. In particular, we show that binomial differential equations without 
linear terms do not have complex centers. Classes of polynomial differential 
equations, with more than two terms, are presented that do not have complex 
centers. We also study the relation between complex centers and the Pugh 
problem. An algorithm is described to solve the Pugh problem for equations 
without complex centers. The method of proof involves phase plane analysis 
of the polar equations and a local study of periodic solutions. 



1. Introduction 
Consider the differential equation 

(1.1) z := J = A N {t)z N + Ak-^z"- 1 + ... + A x {t)z 

where z is complex and Aj(t) are continuous functions in t. Let z(t,c) be the 
solution of (1.1) such that z(0, c) = c. For a fixed real number u>, we say that z(t, c) 
is periodic when z(0, c) = z(tu, c). If the functions Ai(t) are periodic with period u>, 
then a periodic solution of equation (1.1) is a periodic function with period to. The 
multiplicity of a periodic solution ip(t) of (1.1) is the multiplicity of ip(0) as a zero 
of the displacement function q i — > z(uj,c) — c. Note that q is defined and analytic 
in an open set containing the origin. The solution z = is called a center for the 
differential equation if all solutions starting in a neighborhood of are uj— periodic. 
If the coefficients are real functions and Apj(t) does not change sign, then z — is 
not a center (see Q]). There are real centers when A^(t) changes sign (see [3], [5], 
and [15]); the problem is related to the classical center- focus problem of polynomial 
two-dimensional systems. The case N = 3, with A 3 (t) = 1, was considered recently 
in [5J ; where it was shown that the equation could have a center at the origin if the 
coefficients are complex valued. The work in [5J is mainly motivated by a problem 
stated in [5]. It was shown in [5] that if the coefficients are real and bounded by 
a constant C, then there are at most 8exp{(3C + 2) exp[|(2C + 3) n ]} periodic 
solutions. The problem in [S] is to prove an analogue of this result for complex 
equations, that is, for complex z and complex coefficients Ai. This paper is partly 
motivated by the results of [B]. 
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Another motivation for the work in this paper is the local study of the Pugh 
problem about equations with real coefficients. We recall that Pugh problem (see 
|13| ) is to find an upper bound for the number of periodic solutions in terms of 
N. It was shown in [9] that there are no upper bounds for the number of periodic 
solutions for N > 4. Upper bounds can be found for some particular classes; (see 
lj and its references). Hence, one should seek upper bounds for the number of 
periodic solutions in terms of N and the degrees of the real polynomial functions 
Ai(t). A local version of Pugh problem is to find an upper bound for the multiplicity 
of a periodic solution in terms of N and the degrees of the polynomials Ai (t) . This 
problem was considered in [12] , with N = 4. It was conjectured in [12], that if 
A4, = 1, A\ = 0, A2 and A3 are of degree k then the multiplicity of the origin is at 
most k + 3. It was shown in [3], that when k — 2 then the multiplicity of the origin 
is at most 8 and there is a unique equation with this maximum multiplicity. It 
was shown in [2], that the multiplicity of the origin is at most 10 when the degrees 
of A 2 {t) and As(t) are 2 and 3, respectively. Having determined the maximum 
multiplicity, the next step is to construct equations with this number of periodic 
solutions. This is done by making a sequence of perturbations in A2(t) and As(t), 
each of which reduces the multiplicity of the origin by one; a periodic solution 
thus bifurcates out of the origin. This bifurcation task was considered in [5] and 
[3]. The local problem can be considered with the use of Groebner bases method. 
This task is considered in the last section; the case N = 4 and A\ = is studied. 
The problem reduces to study the solvability of system of polynomial equations in 
many variables. Since the solvability in the theory of Groebner bases is over the 
field of complex numbers, it becomes necessary to consider equations with complex 
coefficients. That is to consider complex centers. 

All the known equations with complex centers have linear terms. On the other 
hand, computations for equations without linear part demonstrate that z — is 
not a center for polynomial coefficients. These remarks lead us to conjecture that 
polynomial differential equations without linear terms do not have centers at the 
origin, at least when the coefficients are polynomial functions in t. 

Conjecture. Assume that Aj(t), for j = 2, 3, • • • , N — 1, are polynomial functions. 
The solution z = is not a center for the differential equation 



z 



An^z"- 1 + ■■■ + A 3 (t)z A + A 2 (t)z 2 



First, we present classes of equations with complex coefficients that do not have 
centers at the origin. Our main results in this direction are the following: 

Theorem A. Assume that the differential equation (1.1) has an invariant line 9 = a, 
and either A N (t) > 7 > 0, or A N (t) < 7 < 0. If 7Cos((N - l)a) > then the 
equation can have only a finite number of periodic solutions on the invariant line. 
In particular, z = is not a center. 

Corollary A. Assume that A^(t) > 7 > 0. 

(1) Let Aj(t) = a,j{t) + ibj(t), where aj and bj are real continuous functions. 
If there exists an integer k such that 

2kir 2kir 

GU sm( 7 ) — bi cos( 7 ) 

N — 1 N — 1 

for j = 1, 2, • • • , N — 1, then z — is not a center for (1.1). 
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(2) If the coefficients Aj, j = 1, 2, • • • , N — 1, are real functions then z — is 
not a center for (1.1). 

(3) If N — 1 is a multiple of 4, Aj are pure imaginary functions for all odd j, 
and Aj are real for all even j, then z = is not a center for equation (1.1). 

Theorem B. Consider continuous functions A(t),B(t) and C(t) and let 

Ai = / C(t)dt, 
Jo 

A 2 = / B(t)dt, 
Jo 

and 

A 3 = / [A(t) + (M - L)B{t){ [ C(s)ds)]dt. 
Jo Jo 
The solution z = of the differential equation 

(1.2) z = A(t)z N + B(t)z M + C(t)z L , 

with 1 <L<M<7V = L + M-1, has multiplicity L if and only if A! ^ 0. 
The multiplicity is M if and only if Ai — and A 2 7^ 0; and it is N if and only if 
Ai = A 2 = and A 3 ^ 0. If the origin is a center then Ai = A 2 = A 3 = 0. 

Corollary B. (1) For any continuous functions A(t) and B(t) with A(t)dt 
0, the origin is not a center for the equation 

i = A{t)z N + B{t)z M 

when 1 < M < N. 

(2) Assume that A(t), B(t),C(t), and D(t) are continuous functions, and con- 
sider the differential equation 

z = A(t)z N + B(t)z M + C(t)z L + D{t)z, 

with KL<M<N = L + M-l. Let 

A x {t) = e^- N ^° D ^ ds A{t) 

B 1 {t)=e( 1 - M ^° D ^ ds B{t) 

C 1 (t) = e( 1 -^ti c ^ ds C(t). 

If any of the following conditions is not satisfied then the origin is not a 
center, 

d(t)dt = 0, 



'0 

Bx{t)dt = 0, 

Jo 
and 



/ 

Jo 

f 

Jo 



/ [Ai(t) + (M-L)Bi(t) / C 1 (s)ds]dt = 0. 
Jo Jo 

Remark. It follows from Theorem B, that: 

(1) If B{t) and C(£) are polynomial functions with small coefficients, then z — 

is not a center for the differential equation 

z = z N + B(t)z M + C(t)z L 
with 1 < L < M <N = L + M-l. 
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(2) A polynomial differential equation with only two terms has a center at z = 
only when one of the terms is linear. 

Now, we describe some equations that have centers at the origin. 

Theorem C. The solution z — of the differential equation 

(1.3) z = z N + A(t)z 

is of multiplicity N if and only if 

e J« A{t)dt = 1; f" e (N-l) f* A(s)ds dt ^ Q 

Jo 

The origin is a center if and only if 

e / " A{t)dt = h I"" e (N-l) J* A(s)ds dt = Q 

Jo 

Equations with centers at the origin are given in the following result; the second 
part is a version of the result in [6]. As it was first noticed in [6], the center 
variety in each of the equations considered here contains infinitely many connected 
components. Center variety of each of the known real centers contains a finite 
number of connected components. 

Corollary C. (1) With uj — 2ir, we consider the differential equation 
z = z N + {C + tan(t + co))z, 

where Co is a complex non-real number. The multiplicity of z = is one if 
and only if C = pi, where p is an integer. The multiplicity is N if and only 
if 

/ e( Ar - 1 ^*sec 2 (i + c )d^O. 
Jo 

If this integral vanishes then z — is a center. In particular, there is a 
center when (N — l)p is an odd number. 
(2) Suppose that C is a number and a; is a real number. With lu = 2tt, the 
solution z = is a center for the equation 

i = z N + — — - (C — xi cos t)z 

if and only if C = pi, where p is an integer and a: is a zero of the Bessel's 
function J p . 

Now, we give another class of equations that have a center. The coefficients of 
the equation satisfy a composition type condition similar to those in [4], [5], and 

m- 

Theorem D. Let p(t) be a differentiable oj— periodic function and q(t) is a continuous 
function with e^> I ?( t ) d * = 1. The solution z = is a center for the differential 
equation 

N 

z = q(t) z + £y (t) f k (p(t)) fo z k , 

k=2 

where fk are any continuous functions. 
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Finally, we use the method of groebner bases to study multiplicity of periodic 
solutions when the coefficients are polynomial functions in t, or in cost and sint. 

Theorem E. Consider the following equation 

z = z 4 + A(t)z 3 + B(t)z 2 , 

where A(t) and B{t) are polynomial functions. 

(I) In each of the following classes of coefficients, z — is not a center and the 
maximum possible multiplicity is the same whether the coefficients are complex or 
restricted to be real. 

(1) A and B are polynomial functions in t of degrees 3 and 2, respectively. 

(2) A and B are polynomial functions in f of degrees 5 and 1, respectively. 

(3) A and B are polynomial functions in t of degrees 1 and 3, respectively. 

(II) In each of the following classes of coefficients, z — is not a center and the 
maximum possible multiplicity for real coefficients is less than the maximum mul- 
tiplicity for complex coefficients. Moreover, the origin is not a center. 

(1) A and B are homogeneous polynomial functions in cost and sint of degrees 
1. 

(2) A and B are homogeneous polynomial functions in cost and sint of degrees 
2. 

(3) A and B are homogeneous polynomial functions in cost and sint of degrees 
3 and 1, respectivly. 

In the next section, we prove Theorem A. Section 3 contains the proof of Theorem 
B. The proofs of Theorem C and Theorem D are given in Section 4. In Section 5, 
we consider the case N = 4. Several classes of equations are considered and provide 
evidences to support our conjecture. 

2. Proof of Theorem A 

Proof. Let A m (t) = a m {t) + ib m {t) and z = re 10 . With this notation, equation (1.1) 
becomes 

(f + i9r)e lB = r N e me P N + r^" V^-^ajv-i + i&jv-i) + + ■■■ + re % \a x + ibi). 
Multiplying both sides by e l6 and separating the real and imaginary parts, give: 

JV-2 

r = r N P N cos((A - 1)9)+ ^ r w ^' [ajv-j cos((A- j - l)6)-b N -j sm((N-j - 1)9)} 

3=1 

JV-2 

9 = r^Pw sin((iV-l)0)+ ^^"'[ojv-j wn((N-j-l)6)+b N -j cos((iV-j-l)0)]. 
Now, let us choose a real number a with 



y (N- l) 7 wcos((A- l)a) 
and such that the inequality 

f > Jcos((N-i)a) rN 
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is satisfied when r > a and 9 = a. Consider a solution, re with initial condition 
roe"* . If ro > a, then the inequality implies that 

1 

r(t) > J—. 

(r^~ N - 0.5(N - l)7COs((7V - l)a)t) — 

Therefore, if ro > a, then the solution becomes unbounded at a point to with 

2rl- N 

to ~ (N- l)7oos((JV- l)a) " ^ 

Hence, the solution is undefined on the interval [0, ui]. Therefore, solutions start 
outside the disk r < a are not periodic, and any solution leaves the disk r < a stays 
outside the disk as time increases. 

If there is an infinite sequence of periodic solutions z(t, r n e la ). The numbers r n 
are inside the disk r < a. Hence, let r n tq. If z(t,roe l0! ) is a periodic solution, 
then all solutions starting in a neighborhood of roe la are periodic. But q is an 
analytic function. It follows that q = 0. Therefore, z — is a center for the 
equation (1.2). We define a real number R by 

R = sup{r : z(t, re za ) is periodic}. 

It is clear that R < oo. If the solution z(t,Re la ) is defined, then it is a periodic 
solution. Since q is analytic, it follows from the property of continuous dependence 
of a solution on its initial value, that z(t,re la ) is periodic for r in a neighborhood 
of R. This is contrary to the definition of R. Therefore, z(t, Re la ) is undefined on 
the interval [0,lj] and becomes unbounded at fo£[0, u] and leaves the disk r < a. 
But z(t, Re la ) leaves r < a, whence z(t, re la ) leaves r < a for r < R and close to R. 
These solutions will not return to the disk r < a because f(t) > 0, and hence they 
are not periodic solutions. This is a contradiction to the assumption that z = is 
a center. 

If z(t,roe la ) is undefined, then it leaves the disk r < a. Therefore, z(t,r n e la ) 
leaves r < a for large enough n. This is again a contradiction to the periodicity of 
z(t,r n e la ). □ 

To prove the first part of Corollary A, let a = jrzr- The condition in the 
statement implies that 

dj sin(ja) = bj cos(ja). 
sin((7V -j - l)a) = - sin(a), sin((iV - l)a) = 0. 

Hence, 

9(a) = 0. 

The line 9 — a is an invariant line. Any solution has a point on this line stays on 
the line as long as it is defined. 
On the other hand, 

cos((JV- l)a) = 1. 

Hence, Ajy(t) cos((iV — l)a) = A^(t) > 7. The conditions of Theorem A are 
satisfied and the result follows. 

The second part of Corollary A, follows directly from the first part with k = 0. 
It is known that z — is not a center when the coefficients are real functions (see 
[T]). To prove the third part, we take k = ^j-^- This implies that a — = § ■ 
Thus sin(ma) = if m is even, and cos(ma) = if m is odd. 
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Remark. Equations (1.1), with real coefficients, have been studied in an interesting 
paper of Lloyd [10] using the methods of complex analysis and topological dynamics. 
However, the results of [TU] do not hold when the coefficients are allowed to be 
complex. For example, Abel differential equation, N — 3, has 3 periodic solutions 
when its coefficients are real. On the other hand, it may have an infinite number 
of periodic solutions when the coefficients are complex (see, [6]). However, the 
method of 10J could be modified to obtain global results for equations with complex 
coefficients. This exploitation is more subtle than the real case, and we defer this 
to another paper. 

3. Proof of Theorem B 
Proof. For < t < lj and c in a neighborhood of 0, we write 

oo 

Z(t,c) = J2 d n(t)c U , 
n=l 

where 4(0) = 1 and d n (0) = if n > 1. Thus 

OO 

q(c) = (di(cj) - l)c + J2 d n(u)c n . 

n=2 

The multiplicity is K if and only if 

di(u>) = l,d 2 (uj) = d 3 (uj) =■■■ = d K -i(u) = 0,d K (uj) ^ 0. 

The origin is a center when di(ui) = 1 and d n (uj) = for all n > 1. The functions 
d n (t) are determined by substituting the sum into the equation (1.2) and comparing 
coefficients of powers of c. The following recursive sequence of differential equations 
is obtained with the initial equations. 

d n = ASn + BSm + CSl, 

where 

S K = ^ ' / /: ,/ ,v ' ' • ,/ /.\ • 

jl+32-\ VjK=n 

with Ke{N, M, L}. To solve these equations, we integrate repeatedly. It is clear 
that d\ = and hence d% = 1. The next non-zero equation is 

d\ = C{dx) L . 

It gives 

d L {t) = f C(s)ds. 
Jo 

The formula for du has two possibilities. If M + 1 ^ 2L, then 

d M =B{d l ) M . 

This implies that 

d M (t) = f B(s)ds. 
Jo 

In the case M + 1 = 2L, the equation becomes 

d M ^Bid^ 1 +LCd L {d 1 ) L - 1 ). 
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We integrate this equation to obtain 

d M (t) = J B{s)ds + ^{J^ C{ S )ds) 2 . 

Consequently, the formula for dn has two possibilities. The equation is 

d N = A{d x ) N + BMdLidif 1 - 1 + CLd M (d x ) L - x . 
We integrate this equation to obtain djv(t). If M + 1 ^ 2L, then 

f-t f-t f-S f-t f-t 

d N {t)= I A(s)ds + (M-L) / (B{s) / C{u)du)ds + L / B(a)da / C(s)ds. 
Jo Jo Jo Jo Jo 

In the case M + 1 = 2L, the formula becomes 

d N (t)= f A(s)ds + (M-L) f (B(s) f C{u)du)ds 
Jo Jo Jo 

+L f B(s)ds f C(s)ds+^-{[ C(s)ds) 3 . 
Jo Jo 6 J Q 

The multiplicity is L if f£ C(t)dt ^ 0. The multiplicity is M if f£ C(t)dt = but 
Jq B(t)dt 7^ 0. Finally, if the multiplicity is greater than M, then 

d N (oj)= / A(t)dt + (M - L) I (B(t) / C(s)ds)dt. 



The assumption in the statement of Theorem B, implies that dpf(u)) ^ 0. Therefore, 
the multiplicity is at most N. □ 

The first part of Corollary B follows directly from the above result. If B{t) = 0, 
then c?at(w) > joj ^ 0. Hence, the multiplicity is at most N. 

Now we prove the second part of Corollary B. The first necessary condition for 
a center is e^o D (t) dt = i We make the transformation w = efo D(s)ds z ^ an( j g^ain 

w = A^w" + 5i(t)w M + d(t)w L , 

where A\, B\, and C\ are as defined in the statement of the Corollary. Initial 
conditions and multiplicities of periodic solutions are unchanged under this trans- 
formation. Now the result follows from Theorem B. 

Remark. If the differential equation has a linear term, then the equations for dk are 
more complicated. Instead, we will have linear differential equations in which the 
right-hand side depends also on dk- We consider such cases in the next section. 

4. Proof of Theorems C and D 

Proof. (Theorem C): 

We follow the procedure of the last section. In this case, we obtain a sequence of 
linear differential equations. 

d\ = Adi, 
d\=Ad k , 2<k<{N-l), 

and 

d N = (d!) N + Ad N 
together with the initial conditions 

di(0) = 1; 4(0) = 0, 2 < k < N. 
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Solving these initial value problems, gives 

d 1 (t) = eti A ^ ds , 
d k {t) ee 0, 2 < k < (N- 1), 

and 

d N (t)= [ A(s)ds [ e^-^/o A ^ du ds. 
Jo Jo 

The origin is of multiplicity N if and only if d\{ui) = 1 and ^(cj) ^ 0. These two 
conditions give 

e f£A(t)dt = 1; 

and 



If these two integrals vanish then the solution z = is a center. This follows from 
the general solution of this equation. The general solution is given by 

r t 

w( r 

lo 

where w = z x ~ N and C is a constant. Since w(0) — w(lj), all the solutions are 
periodic. Hence, z — is a center. □ 



,{t) =e (1 - N ^o A ^ ds [(l- N) f^-^o A ^ds + Cl 

Jo 



It is clear that the conditions for a center are not satisfied by any real function 
A(t). We give two equations that have centers at the origin. In the first equation, 
we have 

2 

A{t) = C+ jj—j tan(> + c ). 

Since the zeros of the complex function cos are real, the function A(t) is a periodic 
function of period 27r, when Co is a complex non-real number. The origin is of 
multiplicity one when 

g/o^C+^tantt+co)]* = L 

This gives 

e 27rC eW J r T [ln( S cc 2 (27r+c ))-ln( S cc 2 (co))] _ 1 

Since scc(t + Co) is of period 27r, the condition becomes e 2lTC = 1. Hence, C = pi 
for an integer p. The second condition is 

1-2-K 

e (N-l)pit sec 2( t + c )dt = 0. 

io 

If (N — l)p is odd, then a change of variables 1 1— > t — tt gives 

r 27T 



f" 
Jo 



J 

J TT 



e (N-l)pit gec 2 (t + Co)dt = _ / e (N-l)pU gec 2 (f + Co)dt = Q 

Jo 



Therefore, 

r 2ir 

e (N-i) P it sec 2 (i + co)dt = 0. 

/o 

The conditions for a center given in Theorem C are satisfied. Hence, all the solutions 
z(t, c), with c is in a neighborhood of 0, are 27r-periodic. This proves the first part 
of Corollary C. 



Jo 
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The second part in Corollary C, with N = 3, is similar to the result of [6j. In 
fact, it was shown in [6] that z — is a center for the equation 

z = z 3 + (Co + Cie-** + C 2 e lt )z 



if and only if Co = ki for an integer k, C\Ci ^ 0, and ki\JC\Ci is a zero for 1 fc | - 
Now, we prove the second part of Corollary C. The first condition for a center 
in Theorem C is 

e / 2 " (C-xi sin t)dt _ ^Cir 

But, e 2C7T = 1 if and only if C = pi, for an integer p. The second condition for a 
center becomes 

flit />27T />2tZ 

/ e pit ~ xisint dt = cos(pt- xsint)dt + i sin(pt - x sin t)dt. 
Jo Jo Jo 

The imaginary part of this quantity is zero; it is an integral of an odd 2ir— periodic 

function over the interval [0, 2w\. The real part is 2 J Q cos(pt — xsint)dt. This 

integral is the integral form of the Bessel's function. Therefore The real part is zero 

if a; is a zero of the Bessel's function J p . We recall that the integral form and power 

series expansion of J p are defined by 

J„(x) = / cos(pt — x sin t)dt = > (— lY —. 4 rr- 

p 'Jo f^ jKp + jV- 

The power series representation of Bessel's function is used in [6]. Consequently, 
their proof is much longer than our proof. 

Proof. (Theorem D:) 

With the change of variables z t— > e$o ii s ) ds z ^ the equation becomes 



N 

E 

fc=2 



p'(t)f k (p(t))z k 



From the expansion of Section 3, we can show inductively that the coefficients dk 
are functions of p(t). This implies that z(t, c) is a function of pit). Therefore, in a 
neighborhood of the origin, all solutions are periodic and z — is a center. □ 

Remark. If the coefficient of z N in the statement of Theorem D satisfies the condi- 
tion 

p'{t)f N (p{t)) = e- {N - 1)! ° q{s)ds 
then the coefficient of z N equals 1. This condition is satisfied when p(t) = p 
/jv = 1) and q = 2iri. The period in this case is N — 1. 

5. The Pugh problem 

The use of computer algebra has led to significant progress in the investigation 
of the properties of polynomial differential systems. In this section, we describe 
an application of computer algebra to find the maximum possible multiplicity of 
periodic solutions of polynomial differential equations. In particular, we present 
solutions to the local Pugh problem [T3] and Shahshahani conjecture \¥I . We 
mention that Pugh problem was listed as a part of problem 13 in Steve Smale list 
of 18 open problems for the next century. This problem was considered as a version 
of Hilbert sixteenth problem. Hilbert sixteenth problem is to estimate the number 
of limit cycles of polynomial two-dimensional systems. Research related to Hilbert 
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sixteenth problem has derived enormous benefit from the availability of computer 
algebra as is demonstrated by, for instance, [11] . and [14]. 
Consider the differential equation 

(5.1) i = x 4 + A(t)x 3 + B(t)x 2 

We show that this local problem can be solved using the method of Groebner 
bases. We give an automatic means of finding the maximum possible multiplicity of 
a periodic solution. The algorithm involves computing Groebner bases. Computer 
algebra systems, such as Maple, can be used to implement this algorithm. If the 
equation does not have a complex center, then the local Pugh problem is solvable 
by our procedure. The algorithm for computing the maximum possible multiplicity 
is then described in this section. We apply the algorithm for equations in which the 
coefficients A(t) and B{t) are polynomial functions in t, and in cost and sini. The 
case N = 4 is considered. However, the method works for any N. In the case that 
the coefficients are polynomial functions in t, we show that the maximum possible 
multiplicity is the same whether the coefficients are complex or are restricted to 
be real. When the coefficients are trigonometric polynomials, cases are described 
where the maximums are not equal. Moreover, we show that the origin is not a 
complex center in each of the equations considered; these results provide evidences 
to support our conjecture. 

We follow the same procedure of Section 3. For equation (5.1), d\(t) = 1 and 
the equations satisfied by the d n (t) (for n > 1) are 

(5.2) d n = 2J didjd k di + A 2J didjd k + B 2J 

i+j + k+l=n i+j+k—n i-\-j—n 

These equations were integrated by parts repeatedly and the formulae for d n , with 
n < 8, are given in [3J. The calculations become extremely complicated as n 
increases. It is impossible to accomplish these computations by hand except in the 
simplest cases. 

The next step, in computing the multiplicity, is to consider the quantities 

T] n = d n (u) 

Note thatry„ is a polynomial function in the coefficients of the polynomials A and 
B. The multiplicity of the origin is k if 

r]2 = m = ■ ■ ■ = Wk-i = 0, r\k ^ 0. 

We write [i m ax(C) for the maximum possible multiplicity of z = for equations in a 
class C. For the class of equations in which the coefficients are polynomial functions 
of degree m, Pugh problem is to find \i max in terms of m. 

Now, we give the formulae of rj n with n < 4, for equation 5.1. We use a tilde 
over a function to denote its indefinite integral: 

/(<)= [ f(s)ds 
Jo 

Proposition. For equation (5.1), the quantities 772,773, and 774 are as follows. 

Bdt 

Adt 
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7/4 = CJ + / (BA)dt 

Jo 

Proof. The formulae were obtained by solving equations (5.2) recursively. The 
computations leading to the formulae proceed by sequences of judiciously chosen 
integrations by parts; for any functions / and g, we make use of the identity 

fg = fg- fa- 

These are elementary though together they form a complicated web. To obtain 77„, 
we reduce d n {ui) modulo the ideal generated by d 2 {uj), • • • , d„_i(w). □ 

The Algorithm 

We shall write n max (real) when the coefficients of the polynomials A and B are 
real numbers, and [i max (complex) when the coefficients are complex numbers. It 
follows from the result in [1], that Hmax{feal) < oo. We call the set of equations 
that have this maximum multiplicity, the maximum variety, V ma x- Similarly, we 
define V max (i"eal) and V max (complex). 

To find [i max , first we integrate recursively to compute the functions d n (t) . Then 
we consider the expressions d„(cj), which are polynomial functions in the coeffi- 
cients of A and B. To obtain rjn we reduce d n (uj) modulo the ideal generated by 
d 2 (uj), • • • , d n -\. We stop until the system r\i = 773 = •••=%= has no real solu- 
tions. In this case /i maa; = k. From the theory of Groebner bases, the Groebner basis 
of the ideal (772, 773, • • • , 77^) is (1) if and only if the system 772 = 773 = • • • = 77^. = 
has no complex solutions. So, we have to verify that this maximum multiplicity 
can be attained by certain real values of coefficients. The procedure gives an upper 
bound for \x max . However, for the equations in which the coefficients are polynomial 
functions oit that we will consider, [i max (reaV) — fi max (complex) < 00. 

The algorithm for computing (i max can be summarized as follows: 

• Input: functions A(t) and B(t) which are polynomials in t, or in cost and sini. 

• Integrate to compute d n (t). 

• Compute d n (tu). 

• Find 77„ by reducing d n (uj) modulo (772, 773, • • • , 77„_i). 

• Stop when the Groebner basis (772,773, •• • ,77^.) is (1). 

• Output: \x max — k. For the details related to Groebner bases, we refer to [7]. We 
use Maple8 to compute Groebner bases. 

Now, we consider the first class of coefficients. Here, we assume that u> = 1. 
Let B(t) =d+at + bt 2 and A(t) =C 2 + ct + dt 2 + ei 3 . If 772 = 773 = 0, then 

12Ci + 6a + 46 + 3c = 0, 12C* 2 + 6d + 4e + 3/ = 0. 

This class of equations was used in [5] to construct equations with 10 real periodic 
solutions. We substitute the values of C\ and C 2 from these equations and then 
compute the Groebner basis of the ideal (774, 775, • • • , 77 10 ). This basis is (1). There- 
fore, the maximum possible multiplicity is 10. To find the equations that have this 
maximum multiplicity, we compute the Groebner basis (774, 775, ■ • • ,779). This set of 
equations has 15 solutions, counting multiplicity, and at least one of the solutions 
is real. We summarize the result for this class in the following lemma. 

Lemma 1. If A(t) and B(t) are polynomial functions in t of degrees 3 and 2, re- 
spectively, then ^max(real) — ^, max (complex) — 10. Moreover, V ma x is a zero- 
dimensional ideal. 
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The Groebner basis is computed with the graded reverse lexicographic term order 
and with respect to the list [6, c, d, a, e]; this term order usually gives more compact 
Groebner basis. Then, the basis is changed to the lexicographic term order, which is 
the most suitable to eliminate variables from a set of equations. The Groebner basis 
is given at the end of this section. The last equation in this basis is a polynomial 
in e of degree 15; it has at least one real solution. The other variables are given 
explicitly as functions of e. 

The next class has the coefficients 

B(t) = d+ at, A(t) = C 2 + bt + ct 2 + dt 3 + ei 4 + ft 5 . 

Let 2Ci + a = and 60C 2 + 306 + 20c + 15d + 12e + 10/ = for r] 2 = m = 0. We 
compute the Groebner basis as in the first case. This gives (774, 775, • • • , 7710) = (1). 
Moreover, (774,775,- • • , 779} has a polynomial equation in / of degree 6, which has 
two real solutions and four complex non-real solutions. The other variables are 
given as functions of /. The result for this class are given in the following lemma. 

Lemma 2. If A(t) and B(t) are polynomial functions in t of degrees 5 and 1, re- 
spectively then n m ax(real) = fj, max (complex) = 10. Moreover, V max is a zero- 
dimensional ideal. 

The Groebner basis is given at the end of the section. 
In the last class of polynomials in t, we let 

A(t) = d+bt + ct 2 + dt 3 , B(t) = C 2 + at. 

If r? 2 = 773 = 0, then 

12Ci + 66 + cb + 3d = 0, 2C 2 + a = 0. 

We substitute the values of C\ and C 2 from these equations and then compute 
the Groebner basis of the ideal (774,775, • • • ,77s). This basis is (1). Therefore, the 
maximum possible multiplicity is 8. To find the equations that have this maximum 
multiplicity, we compute the Groebner basis (774,775,- •• ,777). The last equation, 
which is a polynomial in a, has two real solutions and four complex non-real solu- 
tions. The other variables are given in terms of a. Again, the basis is given after 
the statement of Lemma 3. 

Lemma 3. If A(t) and B(t) are polynomial functions in t of degrees 3 and 1, re- 
spectively, then ^J, m ax(real) — fi max {complex) — 8. Moreover, V max is a zero- 
dimensional ideal. 

(8695641600 6 - 6086949120 d - 773773 a 5 + 3151791360 a 2 , 

8695641600c + 13043462400d+ 773773 a 5 - 3151791360 a 2 , 

-50295245 a 4 + 37439568 d 2 + 118513886400 a, 

-120772800 d + 70343 da 3 , 773773 a 6 - 3151791360 a 3 + 3130430976000) 

Now, we consider classes of coefficients which are polynomial functions in cost 
and sin?;. Here, we take w = 2ir. 

Lemma 4- If A(t) and B{t) are homogeneous polynomial functions in cost and sint 
of degree 1 or 2, then fi m ax(real) = 6 and /j, max (complex) = 7. The sets V max (real) 
and V max (complex) are not zero-dimensional ideals. 
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In the case A(t) = ccost + dsini, and B{t) — a cost + bsint. It is clear that 
772 = r)3 = 0. From Maple, we have: 

(?74,f75,f76) = (a 2 + b 2 ,bc -ad- 2). 
This set has complex solutions but does not have a real solution. Moreover 

(m,V5,m, m) = (i>- 

For the other case A(t) = ccos 2 t + dcostsint + c\ sin 2 t and B(t) = acos 2 t + 
b cos t sin t + aisin 2 i. We have 772 = a + a± and ^3 = c + c\. If ai = —a and 
Ci = — c, then Maple gives 

(»74, V5, Ve) = (4a 2 + b 2 , be -ad - 8), 

and 

(V4,V5,V6,V7) = (!>• 

A similar argument proves the Lemma. 

If one of the coefficients A(t) and B(t) contains only terms with of even degrees 
and the other coefficient contains only terms of odd degrees, then it follows from 
the Proposition that fi ma x{real) — fi max {complex) = 4. When the coefficients are 
homogeneous polynomials of degree 3, the Groebncr basis becomes very large. 

Lemma 5. Let A(t) be a homogeneous polynomial of degree 3, and B(t) is a ho- 
mogeneous polynomial of degree 1. The solution z = is not a center. Moreover, 
Hmax(real) = 8, and ^ max (complex) = 10. 

We take B{t) = a cost + bsint. Using the identities sin 3 t = sin£(l — cos 2 1) and 
cos 3 (i) = cosi(l — sin 2 t), we can write a homogeneous polynomial of degree 3 in 
the following form: 

A(t) = c cos t + d sin t + e cos t sin 2 t + f sin t cos 2 t. 

The computer output is given at the end of this section. The Groebner basis of 
(t?4, %, 777) nas a re& l solution; we can take a = d = 1, c = / = 0, e 2 = 432, 
and 36 6 = e. On the other hand, (774, • • • ,fy 8 ) does not have a real solution but 
(774, • • • , 7^ (1), an d (774, • • • , r/io) = (1). On the other hand, ri$ has the form 

m = 7T ((5/8) a 4 + (5/8) 6 4 + (5/4) a 2 b 2 + (1/96) e 2 + (1/96) f). 

It is clear that r\% ^ 0, when the coefficients are restricted to be real numbers. 
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Maple output for Lemma 1. 

(m, ■ ■ ■ , m) = 

(5077631499817345768910374214846690690808088946119335814119700628085731 

205306414251948895671 156736000 a+ 
76568765938297226294215511996032531264938846500344266059142283505 e 14 - 
41366383846391810609354548847965630272631070797484082232995147654360405024 
+9607203693297398146546053698859951537568361672283979690465766480307919 

005421113600 e 8 - 

1077298213530002385272086712469698588685019454221289878562607166668317 

844870406197314322432 e 5 
+5120940699074104547488742022229220044464651558296140078787235311492364 

0199838438141019890846269440 e 2 , 
5077631499817345768910374214846690690808088946119335814119700628085731 

205306414251948895671156736000 6- 
76568765938297226294215511996032531264938846500344266059142283505 e 14 + 
41366383846391810609354548847965630272631070797484082232995147654360405024 
-9607203693297398146546053698859951537568361672283979690465766480307919 

005421113600 e 8 + 

1077298213530002385272086712469698588685019454221289878562607166668317 

844870406197314322432 e 5 - 
5120940699074104547488742022229220044464651558296140078787235311492364 

0199838438141019890846269440 e 2 , 
2542195853515507102698935612213344207862818877038578400901689200003189 

8438696161676309299200 c 
-8114451463043258106360200365311543313307196783429471460595 e 13 + 
4874645150831477756390000057954655294888051734933883104748830146656 e 10 - 
1175665268311675043075984579090225104206445216121635282751902996942202 

131200 e 7 + 

1325591495257443095378482025935592779323248112155531380849978417494633 

00399322103808 e 4 - 
2197311395713698234662544085226573648290742502622945660686244687349190 

8621199707934252072960 e, 
2d + 3e, 

21063213975871712330828163101314609697689436160 e 3 
+3881552635098624688463371219200 e 9 - 16555584308626955352096 e 12 
-438195979208967285885778234437995593728 e 6 + 29908464432145 e 15 
-3407873245453536523174061109544883049290583244800000) 
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Maple output for Lemma 2. 

(»?4, • • • , m) = 



(12699270921091068840199720844433641567807 
091515760 6+ 

2539854184218213768039944168886728313561418303152 e+ 

423247429138101258858306608404799875 f 4 + 
4341516136354618897336295790080428452939652335036/, 
496186409615279857004140496974208217545576225280c- 
595423691538335828404968596369049861054691470336 e- 

38477039012554659896209691673163625 / 4 - 
1028512938658225137299402571203191057331286500736/, 

38477039012554659896209691673163625 / 4 + 
744279614422919785506210745461312326318364337920 d+ 
1488559228845839571012421490922624652636728675840 e+ 
2020885757888784851307683565151607492422438951296 /, - 
76954078025109319792419383346327250 e/ 4 - 
192385195062773299481048458365818125 f+ 
18755846283457578594756510785625070623222781315584000a- 
6849063961008228076202649414049259748983323392 ef- 
17122659902520570190506623535123149372458308480 / 2 , 
4554469328710172509947246108104831527932968138654523985 

40570112 e 2 

+2277234664355086254973623054052415763966484069327261992 

702850560 ef 

28452809860130268863272755268929911307472870565806959073 

22835840 f 

-49539000048633777493317082388741947622188883125 / 5 , 

38477039012554659896209691673163625 f+ 
3984808136948929447103185415887785672672861696 f- 
14218838340357434859762303459719245720157600157874298880000) 



Maple output for Lemma 5. 

(a-l,d-l,c,/,774,--- ,rn) = [-432 + e 2 ,a-l,36&-e,c,d-l,/] 
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(»?4, • • • , m) = 

[fa 4 + 2 b 2 fa 2 + b 4 f -24a 3 + 72 b 2 a, 8 c/a 3 + 16 a 2 fbd + 8 6 3 df + 7 a 2 / 2 6 + 3 6 3 / 2 
-192 a 2 c + 576 abd - 16 6 2 e + 168 afb - 384 6, 24 /a 2 c 2 + 48 d 2 /a 2 + 24 d 2 fb 2 + 42 a 2 df 2 

+ 18 d/ 2 6 2 + 9 fa 2 + 3 / 3 6 2 - 576 ac 2 + 1728 d 2 a - 96 bde 
+4 e 2 a + 1008 adf - 18 fbe + 142 / 2 a - 2304 d - 912 /, 6912 fac 3 + 13824 d 2 /ac + 6912 d 3 fb 
-1536 de 2 6/ + 4 e 3 /a + 12960 dcfa + 9504 d 2 6/ 2 - 406 / 2 6e 2 + 2916 acf 3 
+2676 df 3 b + 31 e/ 3 a + 161 f 4 b - 165888 c 3 + 497664 d 2 c - 82944 ec 2 
+165888 ed 2 - 10368 e 2 c - 384 e 3 + 490752 dcf + 130176 def + 105120 c/ 2 + 23664 e/ 2 , /a 3 

+3 a6 2 / - 24 a 2 + 48 b 2 + 2 6 3 e, -297 ec/ 2 - 297 de 2 / - 1296 dfc 2 - 1296 ed 2 c 
-432 fc 2 - 1296 de/c + 432 ec 3 + 216 e 2 c 2 + 27 e 3 c - 432 e 2 d 2 + 27 df 3 - 52 e 2 / 2 
+e 4 + / 4 + 432 d 3 / + 216 d 2 / 2 , 8 b 2 de + 4 a 2 /c + 12 abdf + 3 b 2 ef + 5 a/ 2 6 - 96 ac + 192 bd 

-8 ae + 64 bf, 96 ed 2 6 - 2 e 3 6 + 48 fac 2 + 144 d 2 /a + 72 debf - e 2 /a + 120 adf 2 
+10 / 2 &e + 23 fa - 1152 c 2 + 2304 d 2 - 192 ce + 16 e 2 + 1536 d/ + 280 f 2 , b 2 e 2 + / 2 6 2 + 12 be 
+12 fa - 288, -2 a/6 + a 2 e - fe 2 e - 24 b, -fb 2 + 2aeb+ fa 2 - 24 a, e 2 a + 12 eac - 12 bde 
-24 adf - 3 /6e - 8 / 2 a - 288 d - 96 /, 12 fac +2e 2 b- 288 c + 24 dea + 9 e/a - 48 e - 12 dfb 

-8 + be - fa - 4 ad + 4 cfe] 

■••,%} = 

[2 e 3 c - 2 de 2 / + 2 ec/ 2 - e 2 / 2 - 2 d/ 3 - / 4 - 11520 a 2 - 23 4 6 2 , 2 e 3 d + 2 /e 2 c 

+/e 3 + 2/ 2 ed 

+2 fc + / 3 e - 11520 ab, e 4 + 2 e 2 / 2 + / 4 + 34560 a 2 + 34560 6 2 , 1440 a 3 
-9 e 2 d - 2 /e 2 - 9 df - 2 / 3 , 1440 a 2 6 + 3 e 2 c + e 3 + 3 c/ 2 + e/ 2 , 1440 6 2 a - 3 e 2 d 
-/e 2 -3 d/ 2 -/ 3 , 1440 6 3 + 9 e 2 c + 2 e 3 + 9 c/ 2 + 2 e/ 2 , -2 a/6 + a 2 e - 6 2 e - 24 6, 
-/o 2 + 2 aeb + /a 2 - 24 a, e 2 a + / 2 a + 288 d + 72 /, e 2 6 + bf - 288 c - 72 e, 
24 ac + 5 ae + 6/, 24 cb + 7 &e + /a - 24, 72 c 2 + 39 ce + 5 e 2 - 3 df - f , 
&e + 24 ad + 7 /a + 24, 24 bd + ae + 5 bf, 24 cd + 5 ed + 5 cf + e/, 72 d 2 - 3 ce + 39 df + 5 / 2 - 
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